M=40; N=1e5;
x_n= randn(1,N);
x_n=2*x_n+16;
figure(1);
hist(x_n,40);
hold on;